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Abstract. The dynamics and the thermodynamics of particles/spins interacting via 
long-range forces display several unusual features with respect to systems with short- 
range interactions. The Hamiltonian Mean Field (HMF) model, a Hamiltonian system 
of N classical inertial spins with infinite-range interactions represents a paradigmatic 
example of this class of systems. The equilibrium properties of the model can be de- 
rived analytically in the canonical ensemble: in particular the model shows a second 
order phase transition from a ferromagnetic to a paramagnetic phase. Strong anomalies 
are observed in the process of relaxation towards equilibrium for a particular class of 
out-of-equilibrium initial conditions. In fact the numerical simulations show the pres- 
ence of quasi-stationary state (QSS), i.e. metastable states which become stable if the 
thermodynamic limit is taken before the infinite time limit. The QSS differ strongly 
from Boltzmann-Gibbs equilibrium states: they exhibit negative specific heat, vanishing 
Lyapunov exponents and weak mixing, non-Gaussian velocity distributions and anoma- 
lous diffusion, slowly-decaying correlations and aging. Such a scenario provides strong 
hints for the possible application of Tsallis generalized thermostatistics. The QSS have 
been recently interpreted as a spin-glass phase of the model. This link indicates another 
promising line of research, which is not alternative to the previous one. 

Key words: Phase transitions; Hamiltonian dynamics; Long-range interaction; Out-of- 
equilibrium statistical mechanics 
PACS: 05.70.Fh,89.75.Fb,64.60.Fr,75.10.Nr 



1 Introduction 

Since the original paper by Ising 1 , magnetic models on a lattice have been extensively used over 
the years to investigate the statistical physics of interacting many-body systems. In particular, many 
generalizations of the Ising model have been proposed. Among them, also models with long-range 
interactions. The thermodynamics and the dynamics of systems of particles interacting with long- 
range forces are particularly interesting because display a series of anomalies with respect to systems 
with short-range interactions '5' . By long-range interaction it is usually intended that the modulus of 
the potential energy decays, at large distance, not faster than the inverse of the distance to the power 
of the spatial dimension. The main reason of the observed anomalies is that systems with long-range 
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forces in general violate extensivity and additivity, two basic properties to derive the thermodynamics 
of a system. Extensivity means that the thermodynamic potentials scale with the system size, i.e. that 
the specific thermodynamic potentials (the thermodynamic potentials per particle) do not diverge in 
the thermodynamic limit. Additivity means that if we divide the system into two macroscopic parts, 
the thermodynamic potentials of the whole system are, in the thermodynamic limit, the sum of those 
of the two components. Though extensivity can be artificially restored by an ad-hoc introduction of a 
N-dependent coupling constant (the so called Kac's prescription 3 ) in the interaction, the problem of 
non-additivity is still present |2|. Thus, the study of long-range magnetic systems on lattices can give 
insights on the statistical properties of long-range potentials, and is therefore useful to investigate the 
statistical mechanics and the dynamics of such systems.^ 

In this paper we study the HMF — a model, a model of planar spins on a d-dimensional lattice 
with couplings that decay as the inverse of the distances between spins raised to the power a 0]|S]. 
When a is smaller than d additivity does not hold and the system shows a series of anomalies. Even 
ensemble equivalence, whose proof is based on the possibility of separating the energy of a subsystem 
from that of the whole, might not be guaranteed in that limit. 

In particular we focus on a = 0. In such a case the model reduces to a mean field model [?ll7ll?ll9l 
I1()II11II12II1H| , since a spin interacts equally with all the others independently of their position on the 
lattice. This case is extremely important since it has been proved that all the cases with a/d < 1 
can be reduced to it. Moreover the dynamical behavior of the model with a = is also a representative 
example of other nonextensive systems 1 4') 1 511 1 611 1 7j . We will discuss the equilibrium phase and 

the dynamical anomalies found in an energy region before the second-order critical point when one 
studies the relaxation towards equilibrium. The connections with Tsallis generalized thermodynamics 
[18ll^l5lllllll9| and with glassy systems pnil21 ' will be also addressed. 

The paper is organized as follows. We introduce the model in section 2, the equilibrium thermody- 
namics is presented in section 3, while the anomalous dynamical behavior and its possible theoretical 
interpretation is discussed in section 4. Conclusions and future prespectives are presented in section 5. 



2 The HMF model 

The HMF — a model has been introduced in ^ and describes a system of classical bidimensional 
spins (XY spins) with mass m — 1. The Hamiltonian is: 

with e = ±1. The N spins are placed at the sites of a generic d-dimensional lattice, and each one 
is represented by the conjugate canonical pair {pi,9i), where the pi's are the momenta and the 9i's 
G [0, 27r) are the angles of rotation on a family of parallel planes, each one defined at each lattice point. 
The interaction between rotators i and j decays as the inverse of their distance rij to the power a > 0. 

Such Hamiltonian is extensive if the thermodynamic limit (TL) A/" — > oo of the canonical partition 
function (In Z) /N exists and is finite. This is assured for each a by the presence of the rescaling factor 
N in front of the double sum of the potential energy. iV is a function of the lattice parameters a, d, N 
which is proportional to the range S of the interaction defined by @E]: 

^^S = y^. (2) 

The sum is independent of the origin i because of periodic conditions. Then, for each a, Va is pro- 
portional to N. When a > d, which we call here the short-range case, S is finite in the TL 5 , and 
things go as if each rotator interacted with a finite number of rotators, those within range S. On the 
contrary when a < d, which we consequently call the long-range case, S diverges in the TL and the 
factor in Q compensates for this. 

^ Very often, as in this case, the term " noncxtcnsivity" is used to refer to the "non-additivity" property. 
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For a = the model reduces to the model introduced originally in Ref.^^ and called HMF 
model. The Hamiltonian of the HMF model is: 

i = l i,j = l 

This model can be seen as classical XF-spins with infinite range couplings, or also as representing 
particles moving on the unit circle. In the latter interpretation the coordinate 6i of particle i is its 
position on the circle and pi its conjugate momentum. For e > 0, particles attract each other or, 
equivalently speaking, spins tend to align (ferromagnetic case), while for e < 0, particles repel each 
other and spins tend to anti-align (antiferromagnetic case) 0. At short distances, we can either think 
that particles cross each other or that they collide elastically since they have the same mass. For 
simplicity of the notation in the following we omit the subscript a = 0. One can introduce the mean 
field vector 

1 ^ 

M = Me"^ = — ^ m, (4) 
1=1 

where rrij = (cos 6i , sin di ) . Here, M and (j) represent the modulus and the phase of the order parameter, 
which specifics the degree of clustering in the particle interpretation, while it is the magnetization for 
the XY spins. Employing this quantity, the potential energy can be rewritten as a sum of single particle 
potentials Vi 

1 ^ 

y = - ^ with v, = l- M cos{0i - (j)) . (5) 

i=l 

It should be noticed that the motion of each particle is coupled to all the others, since the mean-field 
variables M and (j) are determined at each time t by the instantaneous positions of all particles. 



3 Thermodynamics 



In this section we review the equilibrium solution of the HMF model derived in the canonical ensemble 
in ref. [HI . The equilibrium solution of general case HM F — a model can be found in refs. [S] . In the 
canonical ensemble we need to evaluate the partition function: 



Z = j d^pd^e, exp i-PH) 



(6) 



where f3 = l/{kBT), ks is the Boltzmann constant and T is the temperature. The integration domain 
is extended to the whole phase space. It is suitable to factorize the partition function in a kinetic part 



Z 



K 



a Pi exp 



/3 



27r\ 



N/2 



and a potential one 



Zv = exp 



d^6i exp 



-(3e 
'2N 



We shall consider only the ferromagnetic condition, i.e. e = 1. We have 



thus eq. © can be rewritten as 



Z = C I d^e, exp 



(7) 



(8) 



(9) 
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Fig. 1. Temperature T and magnetization M as a function of the energy per particle U in the ferromagnetic case. 
Symbols refer to equilibrium molecular dynamic simulations for N = 10^ and 10^, while the solid lines refer to the 
canonical equilibrium prediction obtained analytically, see text. The vertical dashed line indicates the critical energy 
density located at Uc = 0.75 and f3c = = 2. 



where 





'~I3N' 




2 



(10) 



In order to evaluate this integral, we use the Gaussian identity 

1 



exp 



2 



dyexp[-y2 + y2^x-y] , 
i /i is positive. We can ther« 
d^eJ dyexp[-y2 + y2;iM.y] 

-TT J — OO 



(11) 



where x and y are two-dimensional vectors and /i is positive. We can therefore rewrite Eq. ^ as 

C 



z = - i d'^e. 

IT 



(12) 
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and /i = (3N . We use now definition and exchange the order of the integrals in p2|l . factorizing the 
integration over the coordinates of the particles. Introducing the rescaled variable y yyW/2/3, one 
ends up with the following expression for Z: 



~ o — o / ayexp 



-N{^-\n{2nIoiy)) 



(13) 



where Iq is the modified Bessel function of order and y is the modulus of y. Finally, integral (|13|l can 
be evaluated by employing the saddle point technique in the thermodynamic limit, i.e. for N — s- oo. In 
this limit, the Helmholtz free energy per particle / reads as: 



1- 1^ /27r\ /3 



|^-ln(2^/o(2/)) 



(14) 



(15) 



while the maximum condition leads to the following consistency equation: 

P io{y) 

where Ii is the modified Bessel function of order 1. Eo. ljlSII is the analogous, in the XY model, of 
the Curie- Weiss equation obtained by solving the Ising model in the mean field approximation. For 
/3 < /3c = 2 it presents only the solution y = 0, that is unstable . At (3 = Pc, i.e. below the critical 
temperature Tc = 0.5 (fcs ~ 1), two new stable symmetric solutions appear through a pitchfork 
bifurcation and a discontinuity in the second derivative of the free energy is present, indicating a 
second order phase transition. These results are confirmed by an analysis of the order parameter ^ 

The magnetization M vanishes continuously at Pc (see Fig.^b)). Thus, since M measures the degree of 
clustering of the particles, we have a transition from a clustered phase when /?> /9c to a homogeneous 
phase when fi < Pc- The exponent which characterize the behavior of the magnetization close to the 
critical point is ^ as expected for a mean field model 6 . One can obtain also the energy per particle 
vs temperature and magnetization 

the so called caloric curve, which is reported in Fig. ^a). 

In Fig.l we report temperature and magnetization vs the energy density. The full curves are the exact 
results obtained in the canonical ensemble, the Boltzmann-Gibbs (BG) equilibrium. The numerical 
simulations (symbols) were performed at fixed total energy [microcanonical molecular dynamics) by 
starting the system close to the equilibrium, i.e. with an initial Gaussian distribution of angles and 
velocities, and reproduce the theoretical curves. This is already true for small system sizes as = 100, 
apart from some finite size effects in the homogeneous phase ,2 . As we shall see in the next section the 
scenario is very different when the system is started with strong out-of-equilibrium initial conditions: 
in such a case the dynamics does have difficulties in reaching the BG equilibrium and shows a series 
of anomalies. 



4 Dynamics 

The dynamics of each particle obeys the following pendulum equation of motion: 

N 

e,^^sin{6j-ei)^~Msm{ei-4>) , i = l, ,N (18) 

where M and (j) have a non trivial time dependence, related to the motion of all the other particles in the 
system. Equations (|18|l can be integrated numerically, see for example refs. [SIQ for technical details. 

This is obtained by adding to the Hamiltonian an external field and taking the derivative of the free energy with 
respect to this field, evaluated at zero field. 



6 



Alessandro Pluchino, Vito Latora and Andrea Rapisarda 




t 

Fig. 2. Microcanonical numerical simulations for A'^ = 500 and energy density U = 0.69. In the main part of the figure 
we plot twice the average kinetic energy per particle (which gives the temperature) as a function of time (filled triangles). 
We can easily distinguish a long matastable plateau (QSS regime) preceeding the relaxation towards the Boltzmann- 
Gibbs equilibrium temperature (BG regime). In the BG regime, one finds as expected, a very good agreement with the 
equilibrium thermodynamics value for the temperature, panel (d). In this regime the velocity pdfs reported in panel (b) 
are Gaussians. At variance, in the QSS region, we observe strong deviations from the expected equilibrium temperature. 
Here the specific heat becomes negative, panel(c), and the velocity pdfs, reported in panel (a), are very different from 
the Gaussian equilibrium curve, reported as a full line for comparison. See text for further details. 



In this section we show that, in an energy range from U = 0.5 up to Uc = 0.75, when the system 
is started with strong out-of-equihbrium initial conditions, the model has a non trivial dynamical 
relaxation to the Boltzmann-Gibbs (BG) equilibrium. The class of out-of-equilibrium initial conditions 
we consider, called water hag initial conditions, consists 'm Ot = Q Vi and the momenta uniformly 
distributed (according to the total energy density U). In Fig. 2 we report, for U=0.69 and N=500, the 
time evolution oi 2 < K > /N (where < K > denotes the time averaged kinetic energy), a quantity 
that coincides with the temperature T. As expected, the system does not relax immediately to the 
BG equilibrium, but rapidly reaches a quasi-stationary state (QSS) corresponding to a temperature 
plateau situated below the canonical prediction (dotted curve) . The T vs U plot (caloric curve) , shown 
in inset (c) for the QSS regime, confirms a large disagreement with the equilibrium prediction (inset 
d). Furthermore, it turns out that the system remains trapped in such a state for a time that diverges 
with the size N of the system [H]. This means that, if the thermodynamic limit is performed before the 
infinite-time limit, the QSS become stable and the system never relax to the BG equilibrium, exhibiting 
different equilibrium properties characterized by non-Gaussian velocity distributions (see inset a) [^. 
Such velocity distributions have been fitted in Ref. 9 and have been shown to be in agreement with 
the prediction of the Tsallis' generalized thermodynamics jJHI- The latter is a theoretical formalism 
well suited to describe all those situations where long-range correlations, weak mixing and fractal 
structures in phase space are present jl9| . such as for example excited plasmas ^22^ . turbulent fluids 
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Fig. 3. We plot the time evolution of the variance for the angular displacement, see eq.(18), for A'^ = 500 and various 
energy densities. A ballistic diffusion behavior, 7 = 2 is found for overcritical energies. No diffusion, 7 = 0, is of course 
observed at very small energies, U < 0.25. Anomalous diffusion 7 1.4 is obtained for the case U = 0.6 and U = 0.69. 
This behavior is seen only for a time interval equal to the duration of the QSS regime, indicated by vertical dashed lines 
(see also the previous figure). When relaxation to equilibrium is attained the diffusion becomes again normal, 7 = 1. 
This relaxation time depends on U and diverges linearly with TV 



|28ll24j ■ maps at the edge of chaos j25ll2()ll27ll28j . high energy nuclear reactions ^3 and cosmic rays 
fluxes itOI. 

The characteristics of the QSS have been studied in several papers: vanishing Lyapunov spectrum 
Plll4lll5| . negative microcanonical specific heat^, dynamical correlations in phase-space piTT] . Levy 
walks and anomalous diffusion]?]. Recently also the validity of the zeroth principle of thermodynamics 
has been numerically demonstrated for these metastable states 12 . 

In the following of this paper we focus on the study of anomalous diffusion, long-range correlations 
and the spin-glass phase. 



4-1 Anomalous Dijfusion 

The link between relaxation to the BG equilibrium, and anomalous diffusion was studied in Ref.fTj. In 
that paper the authors studied the variance for the angular displacement defined as 

<^e{t)^jjY.iO.{t)~em? , (19) 

1=1 

which tipically depends on time as u'^(t) oc . The diffusion is anomalous when 7 1 and, in particular, 
it is called subdiffusion if < 7 < 1 and superdiffusion if 1 < 7 < 2. In ref. 7 the authors found the 
existence of an anomalous superdiffusive behavior below the critical energy, connected to the presence 
of the QSS. Superdiffusion turns into normal diffusion after a crossover time Tcrossover that coincides 
with the time Treiax needed for the QSS to relax to the BG equilibrium. We illustrate this behavior 
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in Fig|3| In particular we show for the case N = 500 and various energy densities the time evolution 
for Ug(t). No diffusion, 7 0, is observed for very low energy, while one gets ballistic motion in the 
overcritical energy region, j ~ 2. On the other hand, superdiffusion is found in the energy region 
U = 0.5 — 0.75 with 7 1.4. We plot the cases U — 0.6 and U = 0.69. As expected, however, diffusion 
becomes again normal, 7=1, after complete relaxation. In the figure we indicate with vertical dashed 
lines these relaxation times Treiax for U — 0.6 and U = 0.69. Lines with different slopes are also 
reported to indicate the different diffusion regimes. 

Notice that diffusion starts around U ~ 0.3 when particles start to be evaporated from the main 
cluster j^. 

This result can be interpreted in the following way. When started with water-bag initial condition 
{M = 1) the system immediately decays into the QSS plateau with M w 0. In Fig. 2 this initial part 
is not shown. Then a very slow microscopic relaxation occurs since the force exerted on each single 
spin/particle is almost zero 9 . During the slow relaxation process, many small rotating clusters on 
the unitary circle are continuosly formed and compete each other by trapping the particles in order to 
reach a configuration compatible with the final BG equilibrium state. Thus, particles remains trapped 
for a while in the clusters, until they finally succeed to escape again and so on. Trapping times and 
escape times obey power-law decays with characteristic exponents that allow to relate anomalous 
diffusion with Levy walks. See ref. |7] for more details. Such a mechanism of trappings and Levy walks 
disappears when the system finally reaches equilibrium since at that time only one large cluster is 
present (the system is below the critical point in the ferromagnetic phase). 

Anomalous diffusion can be obtained within a generalized Fokker-Planck equation which generates 
Tsallis distributions with an entropic index q, see refs. |31ll82j . One can extract the following relationship 
between the exponent 7 which characterizes anomalous diffusion and the entropic index q of Tsallis 
thermostatistics 



Considering the value of 7 w 1.38 — 1.4 observed in our case, we would expect then a value of the 
entropic index q w 1.55 — 1.58 which should characterize the dynamical anomalies of the HMF model. 
Actually, this is not the value found in ref. [5] for the velocity pdfs, where we obtained only an effective 
entropic index. However this value is in good agreement to what has been found for the decay of the 
velocity correlation functions which will be examined in the next section. 



4-2 Slow relaxation and aging 

The global effects of the competition between magnetic clusters discussed in the previous subsection 
prevents the system from exploring all the available phase space. Thus one observes a sort of dynamical 
frustration which suggests also an interesting connection with the weak- erg odicity breaking scenario 
typical of glassy systems. Such a scenario, introduced by Bouchaud et al. |2U|, generally occurs when 
the phase-space of the system we consider is not a-priori broken into mutually inaccessible regions, but 
the system remains confined only into a restricted part of it: consequentely one finds slowly decaying 
correlation functions and aging, i.e. the presence of strong memory effects that depend on the history 
of the system. This is just what happens in the HMF model. 

The velocity autocorrelation functions in the QSS regime have been studied in ref|ll). In its simplest 
form the autocorrelation function of the particles velocity can be written as |13| 



1 ^ 

i=l 

If we want to take into account several dynamical realizations (events) in order to obtain better 
averaged quantities, it is possible to use the following alternative definition of the autocorrelation 
function 

r _ <P(t)-P(0)>-<P(t)>-<P(0)> 
crp(rjCTp(Uj 
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Fig. 4. (a) We plot the velocity correlation function Cp(t) vs. time t for the case U = 0.69 and TV = 500 in the 
QSS region, open symbols. The curve plotted as a full line is a q-exponential fit, with q = 1.55. (b) We plot here the 
q-logarithm of the curves and points reported in (a), see text for further details. 



where P = {pi,p2, ■■■Pn) is the N-component velocity vector and the brackets < ... > indicate the 
average over different events, while crp{t) and crp(O) are the standard deviations at time t and at the 
initial time. Such a function has been numerically evaluated in the QSS region. The results are 
shown in Fig. 4 for the case U = 0.69 and N = 1000: a power-law tail is observed, this being a signature 
of long-range correlations. The power-law tail and the initial saturation can be fitted by means of the 
q-exponential relaxation function of the Tsallis' generalized thermodynamics |18|: 



e,(z) = + 



(23) 



with 



Here , t is a characteristic time and q is the entropic exponent which takes into 



account the specific dynamical anomalies of the model into exam. We multiply the curve by a constant 
renormalization factor A (the saturation value). In our case we get for U — 0.69 and N = 1000 a value 
q = 1.55 with a saturation value A = 0.72 and a characteristic time t — 245. By plotting the inverse 
function of the q-exponential, i.e. the q-logarithm defined as 

ln,iz) = ^^-^ , (24) 
1-q 

one can check the quality of the agreement between the theoretical expected behavior and the numerical 
simulations. This is done in Fig. 4(b). 

To study the aging phenomenon we need to evaluate a two points correlation function: 

< Pit + 1^) ■ V{t^) > -< V{t + t^)> ■ < P(i^) > 



(25) 



(Tp{t + tw)ap{tw) 

where t^, is the waiting time. This function has been evaluated in ref. |10llllj . In the QSS region, one 
finds a strong dependence on the waiting time tyj see Fig. 5 (a) for the case U = 0.69 and N = 1000. 
The curves obey a precise scaling as that one observed for glasses. Rescaling the time of the various 
correlation curves by t/t^ with f3 = j, we get a unique curve with a very interesting power-law tail, 
see Fig. 5 (b). Also in this case as in the previous figure, one can reproduce nicely the time evolution 
of the rescaled correlation function with a q-exponential curve. The best fit, reported in the figure, 
was obtained for case q = 1.65, r = 60 and A = 0.7. The q-logarithm is reported in panel (c) in 



10 



Alessandro Pluchino, Vito Latora and Andrea Rapisarda 




ro' 1(f 1(f 10° 70' 1(f 

' t/(tj^'^ 

Fig. 5. (a) For U = 0.69 and N = 1000, we plot the two-times velocity autocorrelation function defined by eq.(25). 
Difli'erent delay times tw are considered, open symbols. Aging, i.e. a marked dependence of the correlation function on tw 
is clearly evident, (b) The curves coUape onto a single one if an opportune scaling is performed, see text. This behavior 
can be reproduced by a q-exponential fit also reported. We plot also in (c) the q-logarithm of the data drawn in (b). 

order to compare the quality of the fit. Notice that this value of q is not very different from the one 
extracted from the decay of the autocorrelation function (22) and from that one extracted from the 
anomalous diffusion. A more detailed investigation in order to deduce analitically the predicted value 
of the entropic index from the observed anomalies is in progress. 

4-3 Spin- glass phase 

In the long-range spin-glass models the aging phenomenon below a transition temperature is associ- 
ated with the complex energy landscape characteristic of the frustrated models. The latter is rich of 
metastable states that play the role of dynamical traps which can confine the system for a long time 
and cause the observed history-dependent slow relaxation dynamics 20 , 3^ . The parallel with spin-glass 
systems can be made stronger by the introduction of a new order parameter for the QSS plateau, the 
polarization p j21j . The polarization represents the temporal average (over a time interval r inside 
the QSS plateau) of the spin vectors, whose modulus is also averaged over the N spins, i.e. 

1 ^ 

i=l 

Such an order parameter allows to characterize in a quantitative way the dynamical freezing and 
frustration of the QSS regime by interpreting the latter as a spin-glass phase. In fact, in analogy 



Dynamics and Thermodynamics of a model with long-range interactions 



11 



with the spin-glass phase of the long-range Sherrington-Kirkpatrick (SK) model |H4j . we have observed 
that, in the QSS regime and in the thermodynamic limit, the magnetization vanishes as M~'^ , but 
the polarization p remains constant around a value equal to w 0.24, see Fig. 2 of ref. (21) . Therefore p 
allows to distinguish between the disorderd QSS glassy phase and the high temperature one. 

This result is very interesting because, at variance with standard glassy systems, neither disorder 
nor frustration are present a-priori in the elementary interactions of the HMF model. On the con- 
trary, they emerge naturally as dynamical features of the model in QSS region and, like the other 
features seen before, disappear when the system reaches finally the BG equilibrium. The introduction 
of the polarization order parameter, which connects the QSS dynamical frustration with the Edwards- 
Anderson order parameter of the SK model, opens a new perspective to understand the true nature 
of metastability in long-range Hamiltonian systems in terms of an emerging glassy behavior. 

5 Conclusions 

In this paper we have briefly reviewed the dynamics and thermodynamics of the Hamiltonian Mean 
Field model. This is a model that allows to study long-range interactions in many-body Hamiltonian 
systems. The equilibrium properties of the HMF model can be derived analitically, though the model 
dynamics presents very interesting anomalies like metastability, superdiffusion, non-Gaussian velocity 
pdfs, long-range correlations, vanishing Lyapunov exponents and weak-ergodicity breaking in an energy 
region below the critical point. These anomalies are common to many long-range systems. Recently 
it has been found that such anomalous regime can be characterized as a glassy phase. Therefore the 
model represents also a very interesting bridge between nonextensive systems and glassy systems and in 
this respect it will deserve more detailed investigations in the future. The link with the nonextensive 
thermostatistics formalism proposed by Constantino Tsallis is also a very promising and intriguing 
line of research. All the anomalies point in that direction and the future years will be crucial in order 
to establish this connection in a firm and rigorous way as it has recently been done for example in 
unimodal maps. 
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